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Abstract. In the paper, we investigate Poincare type inequalities for the 
functions having zero mean value on the whole boundary of a Lipschitz do- 
main or on a measurable part of the boundary. We derive exact and easily 
computable constants for some basic domains (rectangles, cubes, and right tri- 
angles) . In the last section, we derive an a estimate of the difference between 
the exact solutions of two boundary value problems. Constants in Poincare 
type inequalities enter these estimates, which provide guaranteed a posteriori 
error control. 



1. Introduction 

Let f2 be a connected bounded domain in R d with Lipschitz boundary dil. The 
classical Poincare inequality reads 

(1.1) Hkn < Cp(fi)||V™|kn, V^eiJ 1 ^, 

where 

H\n) :={™ 6 1^(0) | lw} a = 0}. 
Here and later on denotes the mean value of g on the set u>. 

It was shown by Steklov [TU] that the constant in is equal to A~ 3 , where A 
is the smallest positive eigenvalue of the problem 

—Aw =Xu in i7; 
( L2 ) d n u=0 on dn, 

It has been shown (see [7]) that for convex domains in R d an upper bound of the 
Poincare constant is expressed throughout the diameter of f2, namely, 

(1.3) C P (Sl) < 

7T 

Other results related to constants in Poincare type inequalities can be found in 
[2 [3l |4j [5j |6] and some other publications cited therein. 

In this paper, we consider estimates similar to for the functions having 

zero mean on a certain part of the boundary (or on the whole boundary) . They are 
as follows: 

(1.4) Hkn<Ci(fi,r)||VH| 2 ,n, Vic e fT^fi.r), 

(1.5) ||tu|| a , r < Ca(n,r)||Vtu|| a ,n, Vu> e H 1 ^, T), 

where T is a measurable part of dtt (we assume that (d— l)-measure of T is positive), 

ff 1 (fi,r) = { l «eff 1 (fi) | Mr = o}. 

Since the quantity ||w||^ := || Vw|j2,o + 1 / r w ds\ is a norm equivalent to the original 
norm of H 1 ^), existence of the constants Ci(Q, T) and C2(ri,r) is easy to prove. 



In this paper, wc find sharp values of the constants in Poincare type inequalities 
for rectangular domains and also for some classes of triangles. Our analysis is based 
on the fact (obtained by standard variational arguments) that the extremal function 
in (|1.4[) is an eigenfunction u G H 1 ^^) of the boundary value problem 

—Ait = Ait in Q; 

^ ^ d n u = \i = 7— / udx on T; d n u = on dQ \T, 

|I I Jn 

which corresponds to the least eigenvalue A > 0. 

Analogously, the extremal function in (|1.5[) is an eigenfunction u € i/ 1 (f2,r) of 
the boundary value problem 

Am=0 in 0; 

(1.7) 

d n u = Xu on T; d n u = on dQ \ T, 

which corresponds to the least positive eigenvalue. 

In both cases the sharp constant in (|1.4[) , (|1.5p is equal to A~ 2 . It is easy to show 
that the eigcnfunctions of the problems (|1.6[) and (|1.7j) form complete orthogonal 
systems in £2^) and in L%(T), respectively. Thus, the analysis is reduced to finding 
the corresponding minimal positive eigenvalues. 

In short, the outline of the paper is as follows. Section 2 is concerned with exact 
constants for rectangular domains in R 2 . In Section 3, we find the constants for 
right triangles and in Section 4 for a parallelepiped. Section 5 is intended to present 
an example, which shows that the estimates can be used in quantitative analysis 
of differential equations. In this example, we consider two elliptic boundary value 
problems with different boundary conditions and source terms. The second prob- 
lem is viewed as a certain simplification of the first one. This means that if the 
functions presenting source terms and Dirichlet or Neumann boundary conditions 
have complicated nonlinear behavior in some sets, then they are replaced by simple 
(e.g., constant) functions. We show that if ft can be decomposed into a collection 
of simple subdomains (for which the constants Cp, C\ and C2 are known), then 
an easily computable bound of the difference between two exact solutions can be 
deduced. We outline that the computation of this bound does not require solving 
a boundary value problem and needs only integration of known functions. In par- 
ticular, this estimate can be used to find a suitable initial mesh in finite element, 
finite difference, or discontinuous Galerkin methods. Our analysis is performed 
with the example of a simple linear elliptic equation. However, by similar argu- 
ments one can obtain similar estimates for other differential equations associated 
with the pair of conjugate operators grad and — div. Other applications of (|1.4[) 
and (|1.5|) are related to a posteriori error estimation methods for partial differential 
equations, where computable bounds between exact solutions and approximations 
often involve constants in Poincare type inequalities (see [8]). 



2. Exact constants for rectangles 

In this section, we assume that Vt is a rectangle with lengthes of sides hi and fo- 
Wc find exact values of the constants in (|1.4j) and ([1.5)1 for the following two cases: 
T coincides with one side of fi and T coincides with the whole <9fl 
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2.1. Case 1: T coincides with one side of the rectangle. In this case it is 
convenient to select the coordinate system such that (see Fig. 1) il = (0, hi)x (0, hq). 
Without a loss of generality, we assume that 

Y = {x x =0,2- 2 e [0,h 2 ]}. 



*2 



£1 

h, 



Figure 1. 



Theorem 2.1. Sharp constants in \l-4\) and fl.5\) are equal to ^ max{2/ii; h 2 } and 
(^tanh(^))"*, respectively. 

Proof. Separating variables we obtain that the eigenfunctions of the problem 
(|F6l) 

/7tto \ /7rfc \ 
u km \x) = cos xi J cos ^-^-^2 j, m = 0, 1, 2, . . . ; k = 1, 2, . . . ; 

/7r(m + i) \ 
"OmW^sm^ - xij, to = 0,1,2,.... 

They form a complete orthogonal system in L 2 (^)- Therefore, the least eigenvalue 

of the problem (| 1 .6[) is min{Aoo; ^10} = mm {(27ir) 2 ' (t^t) 2 }' an< ^ the nrs * statement 
follows. 

Consider another inequality. Similarly, we find that the eigenfunctions of (|1.7[) 
are 

— £2 ) cosh ( — (x\ — hi) ) , k = 0,1,2, 

h 2 / V/12 / 

They form a complete orthogonal system in ^(r). Therefore, the corresponding 

least eigenvalue of the problem (| X . T[) is Ai = jf- tanh(^ il ). □ 

Remark 2.1. It is convenient to present the constant C2 in terms of parameters h 
and n, which characterize the size and the shape oftt, respectively. We set hi = nh 
and ho = h. Then, Co = C*(k)v/i, where C*(k) = , 1 , (see Fiq. 2). 

v ' v ' V 7r tanh 7TK ( y 

2.2. Case 2: T = dil. In this case, the problem is symmetric with respect to 
two axes. Therefore, it is convenient to select the coordinate system such that 
n = (-^,^) x (-If,!f) (see Fig. 3). Due to the biaxial symmetry all the 
eigenfunctions of (|1.6[) and (|1.7|) arc either even or odd with respect to the axes xi 
and X2- 

Theorem 2.2. The sharp constant in is equal to ^ max{/ii; /12}. 
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Proof. First, we consider the eigenfunctions of (|1.6[) . which are odd with respect 
to X\. In this case [i = 0, and we arrive at the following problem: 

, -Au=Xu in n+ : =(0,^)x(-^,^), 

(2.1) 2 ; v 2 ' 2 ; ' 

u=0 on {.Ti=0}nfi, d n u = Q on <9fi+ \ {x x = 0}. 
It is easy to see that the functions 

(i) . /7r(2fc + l) \ (2-Km \ 
U krn\ X ) = 8111 \ ^ X l ) C0S X2 ) ' «,m = 0, 1, . . . ; 

^imW = Sin ( ?r( ' 2 fe i + 1 ' > a:i ) Sin ~ ai2 ) ' fc . m = 0, 1, . . . , 



arc cigcnfunctions to the problem (|2.1[) . They form a system of orthogonal func- 
tions, which is complete in L2(51 + ). Therefore, the least eigenvalue of the problem 

(ED) is A$ = (£) 2 . 

Eigenfunctions of the problem (|1.6[) . which are odd with respect to xi can be 
constructed quite similarly and we find that the corresponding least eigenvalue is 

It remains to consider cigcnfunctions even with respect to both variables. They 
belong to the space # 1 (Q++,r+), where := (0,^) x (0, ?f), and T+ = 

r n <9f2 ++ (see Fig. 3). In this case, we need to solve the problem 

-Au = \u in 

( 2 - 2 ) i j_u A 

d n u = ^ on r+; d n u = on dtt ++ \T+ . 

Moreover, the eigenvalues A£ of the problem (|2.2p (enumerated in the increasing 
order and repeated according to their multiplicity) are critical values of the Rayleigh 
quotient 

«v»ii?„ + * 



(2 - 3) 0W " H|„ + + 

over the space H 1 (n ++ ,T+). 

Consider now the functional Q on the whole space H 1 (Sl ++ ). By the variational 
principle (see, e.g., [T], (1.15)), its critical values A^ enumerated in the increasing 
order and repeated according to their multiplicity satisfy the relation^ A^ < A^ < 
Therefore, if there exists an eigenvalue of the problem (|2.2p in the interval 
(Aq, Af ) then it is necessarily Ag. 

Note that Xf, arc eigenvalues of the conventional Neumann problem 

-Au = Au in d n u = on dfl ++ , 

and thus, Ag = 0, A? = min{(f ) 2 ; (f ) 2 }. 
Now we observe that the equation 

hi /ujho\ h? fuh-\\ 2 

^ Y cot (Y) + f cot (Y) + ^ = 

has a unique solution ujq in the interval (Ojimnl^; j^}) since the function in the 



left-hand side of (|2.4[) decreases from +oo to — oo on this interval. Direct calculation 
shows that the function 

cos(w a;i) , cos(w a;2) 

VniX) = : ; 

V ' sin(^) sin(^) 

solves the problem (|2.2[) with A = Wq. We note that (|2.4p is just the condition 
J r+ v ds = 0. 

Thus, we conclude that Ag = uJq. However, it is easy to see that 

w > mm{ — ; — }. 

hi hi 

Therefore, the least eigenvalue of the problem (|1.6[) is min{(-^-) 2 ; (j^) 2 }, and the 
statement follows. □ 



1 Note that H 1 (n++,r+) has codimension 1 in H 1 (n++). 
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Theorem 2.3. The sharp constant in U.5\) equals (— ^=^=tanh(|£-)) 2 , where 
z = z (a) is a unique root of the equation 

(2.5) tanh ^— ^ tan(zcr) = 1, 



such that zqci < -| , while ag = y "'^j"^ 1 .'^} ■ 

Proof. First, we consider the eigenfunctions of (|1.7|) . which are even with 

respect to both variables. They belong to the space H 1 (fl ++ and solve the 
following problem: 

Ait =0 in 

( 2 -6) 

d n u=\u on r+; d n u = on dfl ++ \ T+. 



Moreover, the eigenvalues of the problem ()2.6|) complemented by zero, enu- 
merated in the increasing order and repeated according to their multiplicity are 
critical values of the Rayleigh quotient 



|V^ll2,n++ 
IMl!,r+ 



over the space H 1 (tt ++ ). Consider another Rayleigh quotient 



5r , _ liv«||» n++ 

Q[v\ 



on the same space. Since Q[v] < Q [i>], by the variational principle its critical val- 
ues A^., which are also enumerated in the increasing order and repeated according 
to their multiplicity, satisfy the relation A| < A|. However, by homogeneity argu- 
ment A| = 2Afc. Therefore, an eigenfunction of (|1.7|) . which is even with respect to 
both variables cannot correspond to the least eigenvalue^. 

Further, we consider the eigenfunctions odd with respect to x\. They lead to 
the following problem in f2 + : 

Au=0 in u = on {xi = 0}; 

(2-7) , 

d n u =Xu on dtt + \ {x\ = 0}. 



We claim that the eigenfunction of (|2.7[) corresponding to the least eigenvalue should 
preserve its sign in Q + . Indeed, the function 

v(xi,x 2 ) = |u(|xi|,x 2 )| -sign(xi) 



2 We can suggest another proof of this fact, which is interesting by itself. Let u be a solution of 
12. 6> . We claim that at least one of sets uj± = C++ D {u 0} has a connected component which 
touches r+ but does not touch the coordinate axes. Indeed, consider a connected component of 
uj-\- touching T+ (in view of the condition J r+ u = 0, such a component exists). If this component 
touches both axes then any connected component of zu- touching T+ is separated either from 
{xi = 0} or from {x2 = 0}. To be definite, let m be a connected component of ro_ which touches 
r+ but does not touch {xi = 0}. Then the function 

v(x 1 ,x 2 ) = u(\x!\, \x 2 \) ■ Xzo(\xi\,\x 2 \) ■ sign(xi) 

belongs to H 1 (f2,T) and provides the same value A of the Rayleigh quotient l|2.8|l as u. If A 
minimizes Q on H 1 (Q, F) then v should be a solution of 11.71 1 which is impossible. Unfortunately, 
this argument is purely 2-dimensional. 
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belongs to 77 1 (ri,r) and provides the same value A of the Raylcigh quotient 



|Vt>" 2 



2.Q 



(2.8) Q[v] = 

ll^ll2,r 

as u. If A minimizes Q on i? 1 (J7,r) then v must be a solution of (jTTTJ) , which 
is possible only if v is positive in + , and the claim follows. Moreover, since 
cigenfunctions of (|2.7[) are orthogonal in L2(dfl + \ {x\ = 0}), an cigenfunction 
positive in f2 + should correspond to the least eigenvalue. 
Now we observe that the equation 

(2.9) tan( — I tanh — 2 - = 1, 



2 / V 2 

obviously has a unique solution cji in the interval (0, |^). Direct calculation shows 
that the function 

v\{x) = sin(wiXi) cosh(wia;2) 

is positive in fl + and solves the problem (|2.7[) with A = lo\ tanh ( " 1 2 / ' 2 ) (the equation 
(|2.9[) is just the equality of quotient d n u/u on sides of rectangle). Substituting zq = 
i Y\/hih2 we conclude that the least eigenvalue of (|2.7p is equal to J^" h tanh(^-) 

where zq is root of (|2.5|) with a = ^/j^- 

In a similar way, considering eigenfunctions of the problem (|1.7p odd with respect 
to a;2 we obtain the least eigenvalue — 2==tanh(=r) where z is root of (|2.5p with 

To complete the proof it suffices to show that the function f(a) = ZQt&nh(^) 
decreases on (0,+oo). We claim that, in fact, af(a) is a decreasing function. 
Indeed, differentiation of (|2.5|) after some transformations yields 

tanh (2a) 

- z oa 



d , ,s 2z (l-tanh 4 (f)) 
(af(a)) 



da K JV " 1 + a 2 -tanh 2 (^)(l-a 2 ) 



1 + tanli (^) 



The fraction here is obviously positive. Further, (|2.5[) implies z a > f. Thus 



tanh(^) 1 7T 

^ ZnQ < < 0, 

l + tanh 2 (f) 2 4 

and the claim follows. □ 

3. Exact constants for an isosceles right triangle 

In this section, we assume that fi is an isosceles right triangle. We find exact 
values of the constants in (|1.4j) and (|1.5[) for the following three cases: T is a leg; T 
coincides with two legs; T is the hypotenuse. 

3.1. Case 1: T is a leg. In this case it is convenient to select the coordinate system 
such that (see Fig. 4a) il — {0 < x 2 < x\ < h} and T = {x\ = h, x 2 G [0, h}}. 

Theorem 3.1. The exact constant in is equal to z^ 1 h, where z is a unique 

root of the equation 

(3.1) zcot(z) + l = 

in the interval (0,7r). 



a 



b 



Figure 4. 

Proof. As in the proof of Theorem [221 the eigenvalues A£ of the problem (|1.6j) 
(which are enumerated in the increasing order and repeated according to their 

multiplicity) are critical values of the Rayleigh quotient Q[v] = „ ir 2 '" over the 

space H^il^). 

Consider the functional Q on the whole space H 1 ^). In accordance with the 
variational principle, the corresponding critical values A£ (enumerated in the in- 
creasing order and repeated according to their multiplicity) satisfy the relation 
A£ < A£ < A£ +1 . Therefore, if the interval (Ap,Af) contains an eigenvalue of the 
problem (|1.6I) . then it is necessarily Aq . 

Note that A£ arc eigenvalues of the conventional Neumann problem 

(3.2) — Am = Xu in f2; d n u = on <9fL 

By even reflection with respect to the line {x\ = X2] we conclude that any eigen- 
function of (|3.2[) is an cigenfunction of the Neumann problem in the square (0, h) x 
(0,h). In particular, Aq = corresponds to the eigenfunction uq = 1, and Af = 
(^) corresponds to the eigenfunction u\(x) = cos( 2 ^ 1 ) + cos(^|p-). 

Now we observe that the equation (|3 . 1 [) obviously has a unique solution in the 
interval (0,7r). Direct calculation shows that the function 




solves the problem (|1.6p with A = (4^ ) 2 (the equation (|3.1|1 is just the condition 
J r vods = 0). Thus, we conclude that Aq = (^;) 2 , and the statement follows. □ 

Remark 3.1. Approximate value of the root in k3.1\) is 2.02876. Thus, the constant 
in Theorem \3.1\ is approximately 0.4929/i. 

Theorem 3.2. The sharp constant in U.5\) is equal to tanh(zo)) 2 where Zo 
is a unique root of the equation 

(3.3) tan(z) + tanh(z) = 

in the interval (0,7r). 

Proof. Using the monotone rearrangement (see, e.g., [5]) with respect to X2 
we can suppose that the minimizer v of the Rayleigh quotient Q[v] over the space 
JT 1 r) is monotone decreasing in X2- Therefore, w|r has exactly one change of 
sign. Moreover, since any other eigenfunction of the problem (| 1 . T[) is orthogonal to 
1 and to v in ^(r), an eigenfunction u such that w|r has exactly one change of 
sign should coincide with v up to a constant multiplier. 



Now we observe that the equation (|3.3[) obviously has a unique solution in the 
interval (0,7r). Direct calculation shows that the function 



vi(x) = cos y— — J 



i z x 2 \ , , fz xi\ fz x 2 
cosh — ; — + cosh — ; — cos 



h 

solves the problem (|1.7[) with A = ^ tanh(zo) (the equation (|3.3[) is just the condi- 
tion J r vi ds = 0). Since ^ < 2b < n, «i|r is monotone decreasing in x 2 . Therefore, 
Vi\y has exactly one change of sign, and the statement follows. □ 



Remark 3.2. Approximate value of the root in h3. 3\) is 2.3650. Thus, the constant 
in Theorem \3.2\ is approximately 0.6560/i2. 



3.2. Case 2: T coincides with two legs. In this case we again assume that (see 
Fig. 4a) fl — {0 < x 2 < xi < h} and r = {x\ = h, x 2 € [0, h]}U{x 2 = 0, x\ € [0, h}}. 

Theorem 3.3. The sharp constant in {l-4ty is equal to 

Proof. By even reflection with respect to the line {x\ = x 2 } we conclude that 
any eigenfunction of (| 1 .€>[) is an eigenfunction of the same problem in the square 
fi' = (0, h) x (0, h) with T = 9fi'. This problem is solved in Theorem [2~2l and the 
least positive eigenvalue equals (^) 2 - The dimension of corresponding eigenspace 
equals 2 and contains the function cos(^xi) + cos(^x 2 ) which solves the original 
problem in the triangle. □ 

Theorem 3.4. The sharp constant in U.5)) equals (22a tanh(z )) 2 , where zq is 
a unique root of the equation \2. 5\) with a — 1 such that z < ^ . 

Proof. We again use even reflection with respect to the line {x\ = x 2 } and 
reduce our problem to the problem in the square (0, h) x (0, h). This problem is 
solved in Theorem 12.31 and the least positive eigenvalue equals 2£a tanh(zo)- The 
dimension of corresponding eigenspace equals 2 and contains the function 

sin ( Zo (^i _ i)) cosh ( Zo (^e _ i)) + cosh (, (^i _ i)) sin ( Zo (^a _ i)) 

which solves the original problem in the triangle. □ 

Remark 3.3. Approximate value of the root in A2.5\) with a = 1 is 0.93755. Thus, 
the constant in Theorem \3.2\ is approximately 0.8523ft.2. 

3.3. Case 3: T is the hypotenuse. In this case it is convenient to select the 
coordinate system such that (see Fig. 4b) Q = {0 < \x 2 \ < Xi < h} and T = {x\ = 
h,x 2 € [—h, h]}. 

Theorem 3.5. The sharp constant in is equal to Zq h, where Zq is defined 

in Theorem \3.1\ 

Proof. First, we consider the eigenfunctions of (jl.6|) . which are odd with respect 
to x 2 . In this case \i = 0, and we arrive at the following problem in fi + = {0 < 
x 2 < x\ < h}: 

—Am =Xu in f2 + ; u = on jx? = 01; 
(3.4) J 12/, 

d n u=0 on dil + \ {x 2 = 0}. 

Similarly to Theorems 13.31 and 13.41 we use even reflection with respect to the line 
{xi = x 2 } and reduce (|3.4j) to the problem in the square (0,h) x (0,h). Thus, 
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wc conclude that the least eigenvalue of the problem (|3.4j) is equal to \{j L ) 2 and 
corresponds to the eigenfunction 



u (x) = sm [—) sm [— 



Next, we consider the eigenfunctions, which are even with respect to x 2 . Then we 
arrive at the problem (| 1 .€>[) in f2 + which is solved in Theorem 13. II 

To complete the proof we compare 5(f) 2 and (jP- ) 2 . It is easy to check that 
■ cot(-^=) < — 1. Since 1 i-> t ■ cot(i) is a decreasing function on (0, tt), this means 
-7= > Zo, and the statement follows. □ 

Theorem 3.6. The exact constant in hl.5\) is equal to hi. 

Proof. First, we consider eigenfunctions of the problem (|1.7[) even with respect 
to x 2 . Then we arrive at the problem (|1.7[) in Q + which is solved in Theorem 13.21 

Further, let us consider the eigenfunctions, which are odd with respect to x\. 
We arrive at the following problem in 

Au = in d n u = on {x\ = x 2 }; 

(3.5) 

u = on {x2 = 0}; d n u = Xu on {x\ = h}. 

Direct calculation shows that the function x\x 2 is positive in fi + and solves the 
problem (|3.5p with A = j- . Similarly to the problem (|2.7j) , it should correspond to 
the least eigenvalue. 

To complete the proof we compare # tanh(zo) and \. Since 2"o > 5, we have 
zotanh(zo) > §tanh(^) > 1, and the statement follows. □ 

4. Constants in three dimensional case 



Theorems I2.1H2.3I can be extended to functions of three variables. The corre- 
sponding proofs are quite similar. Therefore, we present them in a concise form 
paying major attention to 3D specifics. 

Theorem 4.1. Let ft = (0,/ii) x (0,/i 2 ) x (0,h 3 ) and T = <9ft n {x x = 0}. 
Then the sharp constants in {l.J$ and U.5\) are equal to — max{2/ii; h 2 \ h 3 } and 

V max fen tanh ( ma x{fc;M )) 2 ' respectively. 

We omit the proof, which is quite similar to the proof of Theorem 12. II 

Remark 4.1. Let h = ma,x{h 2 ; ^3} and hi = nh. Then Ci(ft, T) = — max{l, 2k} 
and C 2 (0,r) = C*(k)V^, where C*(k) = , 1 = (see Fig. 2). 

Theorem 4.2. Let ft = (-^-, \) x (-4f , \~) x (-\, ?f) and T = <9ft. Then the 
exact constant in {l.J$ is equal to ^ max{/ii; h 2 ; h 3 }. 

Proof. The proof is similar to the proof of Theorem 12.21 Instead of (|2.4p we 
obtain the equation 

h\h 2 / wh 3 \ h-\h 3 /uih 2 \ h 2 h 3 /coh^\ 2 



■ cot 



Its unique solution in the interval (0, (here h = max{/ii; h 2 ; h 3 }) is greater than 
x> and the statement follows. □ 
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Theorem 4.3. Let and T be as in Theorem Assume (for the sake of def- 
initeness only) that hi < h% < h%. Then the exact constant in il.5\) is equal to 

tanh(zi)) 2 , where (z\,Z2) is a unique solution of the system 

Z\ z% 
(4.1) — tanh(zi) = — tanh^) = 
h\ h-2 




tanh (zi) ( Zih% / tanh (z\ 

■ cot — 4/1 



tanh 2 (z 2 ) V hi V tanh 2 (z 2 ) 



such that ^0 



tanh 2 (zi) ,r 
tanh a (z 2 ) 2 ■ 



Proof. Similarly to the proof of Theorem 12.31 we conclude that the eigenfunc- 
tion of (|1.7[) corresponding to the least eigenvalue must be odd with respect to one 
of the coordinate axes ( for instance, with respect to X3). This gives the following 
analog of (377) in the domain h+ = (-&-, !f) x (-ff, !f) x (0, 

Au =0 in u = on |x 3 = 0|; 

(4.2) J 1 " 

9 n u =Au on c?Sl + \ {0:3 = 0}. 

Repeating the proof of Theorem 12.31 we find that the eigenfunction of (|4.2[) , which 
is positive in Sl + , corresponds to the least eigenvalue. 
It is easy to see that the equation 

( lihi\ fvh-2 



(4.3a) /itanh y— — J = i/tanh ^ 

defines an increasing function v = v{jj) on R + , and the equation 

(4.3b) /itanh (^i) = vV + ^) ■ cot (v^T^O y) 



has a unique solution /io such that yMo + 1/2 (Mo) 2 < f ■ Direct calculation shows 
that the function 



?7i(ar) = cosh(/i xi) cosh(^(^ )^2) sin f y / /x^ + i/ 2 (//o)x3 



is positive in f2 + and solves the problem (|4.2p with A = /^otanh (^V 1 -) (note that 
(|4.3[) reflects the boundary conditions on the sides of parallelepiped). Substituting 



zi = z 2 = J/ ^ t °' fl2 ; we conclude after some manipulations that the least 

eigenvalue of (|4.2j) equals 2 ^ L tanh(^i) where (21, £2) is solution of (|4.1[) . 

The eigenfunctions odd with respect to other coordinates can be constructed 
quite analogously. However, some additional calculations show that Ui is the best 
eigenfunction provided that /13 is the longest edge of f2. Thus, the statement follows. 
□ 



5. An application of the estimates 

In this section, we discuss the meaning of the above derived estimates for quan- 
titative analysis of solutions to partial differential equations with the paradigm of 
a linear elliptic equation. However, similar analysis can be performed for other dif- 
ferential equations associated with the pair of conjugate operators grad and — div. 
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Assume that the boundary dfl consists of two measurable nonintersecting parts 
T D and T N associated with Dirichlet and Neumann boundary conditions, respec- 
tively. Consider the following elliptic boundary value problem V: 



(5.1) div(^Vw) + / = in ft, 

(5.2) u = u on T D , 

(5.3) AVu-n = F on T N . 
Here the dot stands for the scalar product of vectors, 

(5.4) feL 2 (n), FeL 2 {T N ), and uq e H 1 ^). 



We assume that the matrix A is symmetric, bounded, and satisfies the uniform 
cllipticity condition 

M-Z > c|?| 2 , c>o, 

Standard (generalized) solution to the problem V is a function u £ _ff 1 (fi) such 
that 

u — Uq £ Vo := {w £ H 1 ^) : w\ r o = 0} 

and 

(5.5) J A\7u-\7wdx = J fwdx+ J Fwds, £ V . 

Well known results in the elliptic theory guarantee the existence and uniqueness of 
the solution it. 

In practice, finding a solution u is replaced by finding a sequence of approximate 
solutions u n converging to u. Usually, u n is constructed as a Galerkin approxima- 
tion associated with a certain finite dimensional space V n If the functions /, F, 
and Uo are complicated (e.g., rapidly change or oscillate in some parts of CI, T N , 
and T D ), then finding u n may be a difficult problem. For example, if approxima- 
tions are constructed with the help of simple (e.g., piecewise affine) functions and 
the boundary conditions are defined by complicated nonlinear functions, then the 
boundary condition on T D and T N cannot be exactly satisfied. Similar difficulty 
arises if a curvilinear boundary is approximated by piecewise affine functions. Nu- 
merical computation of the integrals involving / and F (which is necessary in all 
variational-difference numerical schemes) leads to errors generated by the fact that 
on mesh cells the source terms are usually simplified (e.g., replaced by mean values). 
All these errors have a common source: they arise because in reality the construc- 
tion of u n is based on a different problem P. Approximation and integration errors 
induce additional errors in discrete solutions, which are usually estimated only in 
an asymptotic sense. However, in quantitative analysis we need concrete values of 
them. Indeed, if such an error is smaller than the desired accuracy level (which 
in the majority of cases is known a priori), then we can ignore inconveniences in 
boundary conditions and inaccuracy in local representations of source terms. On 
the other hand, if it is essentially larger, then the mesh and integration methods 
are invalid for our purposes (i.e., this is a signal to use a finer mesh and/or more 
accurate integration methods). Thus, a guaranteed and easily computable estimate 
of the error can help to select suitable meshes, approximations of source terms, and 
quadrature formulas without directly solving a boundary value problem. 

Below we show how the required estimate can be deduced with the help of the 
Poincare type inequalities considered in previous sections. We note that estimates 
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of errors caused by simplification of the coefficients entering A has been recently 
derived in [§] , so that estimation of summed effect can be done by combining these 
estimates. 

Let us assume that f2 is split into a set O of " simple" nonoverlapping subdomains 
fli (e.g., they can be cells of a certain mesh). Each J7j belongs to one of the following 
three subsets: 

o D -.= {fljcn | on, n r D ■= rf ^ 0}, 

O n := {QjCQ | dfli n T N := rf + 0}, 
O 1 := 0\(0 D UO N ). 

In other words, O 1 contain interior subdomains, O d contain subdomains associated 
with F D , and elements of O n are the subdomains associated with T N . Then, 
O = H D U H 7 U Ti N , where fl D ,n N , and fl 1 consist of Q< from O d , O n , and O 7 , 
respectively (see Fig. 5). 
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Figure 5. 



Now, instead of V we consider a modified (simplified) problem V: 

(5.6) div(lVu) + /=0 in Q, 

(5.7) u = u on T D , 

(5.8) AVu-n = F on T N . 

Our goal is to deduce an estimate of the difference between the exact solutions of 
these two problems. For this purpose, we define the energy norm 

\\u — u||| 2 := J AV(u — u) ■ V(u — u) dx 
n 

and note that, by the ellipticity condition, 

(5.9) > Vc- \\V(u-u)\\ 2 ,n. 
Now we introduce the quantities 

1/2 

(5-10) Dx~ ( J2 C?||/-/lll in< ) , 

n,eo 

(5.H) D 2 := (J2 C^TffWF-FWl^) 112 , 
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(5.12) 



where 

c x {n h TP) if si,eo fl 

These quantities are easily computable provided that the constants Cp, Ci, and C2 
associated with the corresponding subdomains are known. Indeed, if /, /, F, and F 
arc defined, then finding the quantities is reduced to integration of known functions 
and does not require solving a boundary value problem. The theorem below shows 
that a guaranteed and directly computable bound of the difference between two 
exact solutions can be expressed throughout D\, D 2 and other easily computable 
quantities. 

Theorem 5.1. Let u and u be the solutions of i5.1)) - ^573\) and i5.6]) - [57E\) , re- 
spectively. Suppose also that the right-hand sides of i f 5. 6\) and i5. 8\) satisfy the 
conditions 

(5.13) {{/-/Sf2,=o vaeo'uo", 

(5.14) {{F-FS r «=0 VSJ.eO" 
Then 



(5.15) |u - u|| < pi + ^p 2 + Pi 
where 

(5.16) 2 Pl = Dl + ° 2 + 101; p 2 = 1 +l 1 {cf>) +I 2 (<f>), 

(5.17) 2o= E lu -u } r n f(f-f)dx, 

(5.18) l 1 (<f>)=f(f-f)cj>dx, M4>)= f(F-P)4>ds, 

n r N 

and <fi is an arbitrary function in H 1 ^) such that <fi = uq — u~o on T D . 
Proof. We use (I5~5l). (I5~6l)-(l5~8l) and obtain 



(5.19) J A(Vu - Vu) • Vw dx = f{f-f)wdx+ J(F-F)wds Ww G V . 
a n r N 

Since w = u — u — G Vo we can use it as a trial function. Then 

(5.20) J AV(u~u) -Vwdx =\\u-u\f + J AV(u - u) ■ V<j)dx > 



> Ml — u|| 2 



Consider the first term in the right-hand side of (15. 19ft : 
(5.21) f(f-f)wdx 



E f(f-f)(u-u)dx+ J2 /(/-/)( 



u — u)dx +Ii((f>). 
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The terms of the first sum in (|5.21[) arc estimated using (|5.13l) and (|f .f I) : 

(/ - f)( u - u)dx = (f - f)(u -u- fu-ujajdx 



<C , p(n i )||/-/l| 3 ,n t ||V(«-u)|| 2 ,n i , n i eO I UO N . 
The terms of the second sum are estimated with the help of (|1.4p as follows: 

(/ - /)(" ~u)dx= / (/ - f)(u -u-{_u- u} r o)dx 



Si, 



+ / (/ - /){{ « - uJ r o dx < Tf )||/ - /|| 2 , n , ||V(« - u)| 



2,f2i 



+ {{u-u} r n (f-f)dx, OiGO 



Summing up these estimates and using (|5.9|) and (|5.17|) we obtain 
(5.22) f[f-f)wdx< Y, C 4f~ fh,n t \\ V(«-u)||a inj +lo+M<f>) 



n t eo 



In a similar way, using (|5.14|) and (|1.5j) we deduce 

(5.23) J(F - F){u -u)ds< C 2 (fi, rf )||F - F|| 2 ^ || V(u - u)||a,n«, 

and, therefore, by (|5.9[) and (|5.18p . 

(5.24) y (F-F)wds< %u-0||+J 2 (# 

Now ([S"Tl?j) , flOD, dOg) , and ([5~2^|) imply the estimate 

(5.25) |w-0| 2 < 2pi|«-u|| +p 2 , 

where the quantities p\ and p 2 arc defined by (|5.16[) . 

The quadratic inequality (|5.25|) easily implies (|5.15l) . □ 

Theorem 15.11 presents the most general form of the estimate. If uq = Uo, then 
this estimate can be significantly simplified. Indeed, in this case one can choose 
0=0, and (|5 . 1 5[) is reduced to 

(5.26) ||| U _£|||<___L__. 

Moreover, in this case we can replace C\(£li,rf) in (|5.12|) by a smaller constant 
C F (Cli,TP) such that 

Hkn < C P (Sli,T?)\\Vwh,n, Vw e H 1 ^) : w\ r o = 0. 

For simple domains such as rectangles or isoceles right triangles this constant is 
well-known. 

is 
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